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ABSTRACT 


The evolution and long-time stability of a double layer (DL) in a discrete auroral arc requires that the parallel 
current in the arc, which may be considered uniform at the source, be diverted within the arc to charge the flanks of 
the U-shaped double layer potential structure. A simple model is presented in which this current redistribution is 
effected by anomalous transport based on electrostatic lower hybrid waves driven by the flank structure itself. This 
process provides the limiting constraint on the double layer potential. The flank charging may be represented as that 
of a nonlinear transmission line. A simplified model circuit, in which the transmission line is represented by a 
nonlinear impedance in parallel with a variable resistor, is incorporated in a one-dimensional simulation model to 
give the current density at the DL boundaries. Results are presented for the scaling of the DL potential as a function 
of the width of the arc and the saturation efficiency of the lower hybrid instability mechanism. 


I. INTRODUCTION 


A vast body of ground-based, rocket, and satellite observations reveals that auroral-zone acceleration 
processes occur in a hierarchy of latitudinal scale widths. On the scale of the inverted-V region (AA > 1°) parallel 
electric fields are observed in narrow, soliton-like structures interpreted as weak or ion-acoustic double layers 
(DL’s) (Temerin et al., 1982). Assuming statistical homogeneity of the distribution of these weak DL’s over an 
altitude range comparable to 1 Re, one infers a total potential drop of up to a few kV, typical of the inverted-V 
region. On smaller spatial scales (AA <0.1°) more energetic precipitation is observed in discrete arcs, which have 
projected widths ~1 km in the ionosphere. Discrete arcs (DA’s) are associated with electrostatic shocks (Torbert 
and Mozer, 1978; Kletzing et al., 1983). We adopt the hypothesis that electrostatic shocks constitute the nearly 
field-aligned “flanks” of the paradigmatic U-shaped potential structure of a strong double layer. Although this 
hypothesis seems plausible, many questions exist concerning the conditions under which DL’s may exist in space, 
their dynamics, and their structure. These questions are vital for understanding the complex observational morphol- 
ogy of fields and particles in the auroral zone. At present, investigations of such questions must to a large extent be 
motivated by and proceed from consideration of the fast-growing literature on experiments and simulations, 
although usually the applicability of these situations to DL’s in space is indirect (Smith, 1985, 1986a). 

In this paper, we first discuss theoretically the question of what limits the potential of DL’s in auroral arcs, 
and report results of recent simulations of DL’s in a model circuit. Somewhat more detailed expositions are given by 
Smith (1986b, c). 


PRECEDING PAGE BLANK NOT FILMED 


211 


II. THEORY 


Experiments and simulations (Goertz and Joyce, 1975; Coakley and Hershkowitz, 1979) reveal a scaling 
law for the DL potential in terms of its length fi DL and the electron density n^ on the low potential (cathode:K) side; 
we write this law as 


where £ D l = <j>Diy m ax|E|||. In space, £ DL is not limited a priori and, absent other constraints, equation (1) implies 
that the potential may grow to much larger values than the observed limit on the auroral precipitation energy, which 
is a few tens of keV. 

This dilemma is resolved by considering how the field-aligned flanks of the arc become charged during the 
evolution of the DL’s. We adopt as a starting point the basic idea of the recent MHD models discussed by Haerendel 
( 1983 ) and Goertz ( 1985 , 1986 ) in which the DL evolves in the parallel current sheet of a kinetic Alfven wave. This 
scenario limits the thickness of the sheet a priori to a few times the ion gyroradius at an energy representative of the 
distant plasma population in the generator region. Taking this energy to be — 1 keV, we may estimate j| t by 


in = neU - ~ ^ 


8B i 

rT 


Using n ~ 10 2 cm 3 , B 0 ~ 0.05 G, SBj^ ~ 10' 3 G, and assuming a current sheet of a few kilometers thick, this 
equation gives a relative drift velocity U ei greater than the electron thermal velocity V e . At such a relative drift 
velocity, the current sheet is unstable to a variety of instabilities, including the ion cyclotron and Buneman 
instabilities. We expect the instability to be triggered at some altitude z* where the density and magnetic field 
profiles first combine such that U ei exceeds the threshold drift. In addition, experiments reveal that the U-shaped 
structure, with the field-aligned flanks curved toward the low potential side as is required for Earthward-directed 
Poynting flux (Smith, 1 986a), requires o> e < fl e , where o> c and fl e are the electron plasma frequency and gyro- 
frequency, respectively. This is just the condition for strong magnetization (io e /D e = R e /\ e ), and is fulfilled in a 
limited altitude range along the auroral field lines (Gumett, 1974). 

Simulations show that DL’s evolve from current-driven instabilities when the current is interrupted by trap- 
ping (Smith, 1 982a, b). Trapping creates local regions of macroscopic non-neutrality; in the finite-thickness current 
sheet, the plasma tends to expel charge in the transverse direction in an attempt to neutralize the local electric field 
(Fig. 1). Electrons are tightly magnetized and cannot be expelled very far, but the ion motion is essentially ballistic 
(the evolution time scale is <fV), and ions are accelerated in the transverse direction out to some distance greater 
than their gyroradius. Owing to mirror forces, the expelled charge spreads upward, providing the initial charging of 
the flanks. 

The charging mechanism described above operates in the transient phase. The characteristic time scale of the 
evolution is t dl ~ ^Di7 UjA> where Uj A > C s the ion inflow velocity in the frame of the DL. The charge spreads along 
the flank at velocity c/A/e^ , where £j is the dielectric constant. In the MHD limit, tj — c 2 /V A 2 , but we shall see later 

that € j ‘ c ; v A in tuc ul uank. In any case, however, we find that the time Tj nd = /dsc/V ,/z /c for the charge to 

spread along the field lines back to the generator region (Fig. 2) is long compared to ioi, < t < T jnd , the spreading 
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charge would thin out along B and the flank would not sustain the DL potential; then the DL would discharge. This 
would occur in a time short compared with the typical lifetime of discrete arcs. Therefore, 
asymptotic stability of the DA requires a transport mechanism producing a cross-field current density J x (4>), which 
persists as d<J> DL /3t -*-0. 

If such a mechanism exists, then in the time-asymptotic regime, the density n cK in equation ( 1 ) is determined 
by current continuity and is of the form 


^eic(z*) f*eoo § J"z. dz J X (4>,Z) , (2) 


where g is a constant. Substituted in equation (1), equation (2) provides the physical constraint on the DL potential, 
which is transparently self-stabilizing for J x , a monotonically increasing function of <j>. 

A mechanism to maintain a distributed J x in the time-asymptotic regime is discussed by Smith ( 1 986b) . The 
mechanism is based on anomalous transport due to lower hybrid waves which are driven by the inhomogeneous 
structure of the flank itself. The discussion above implies that the initial scale length of the perpendicular electric 
field E x in the flank is 2f > Rj. As this field is established along B, the electrons acquire the local polarization drift 
velocity cE x /B 0 . The ions, however, encounter an inhomogeneous electric field over the scale of their gyroradius, 
and so their drift orbit is modified by finite-Larmor-radius (FLR) effects. For 8 f > R i; the ion drift speed is approxi- 
mately given by the first non-vanishing order of the phase-averaged FLR correction: 

>(>-#)s$ ,3) 


Then there is a relative drift 



if U ej > Vj, this relative drift drives the electrostatic modified two-stream instability (MTSI) studied by McBride et 
al. (1972). (Other instabilities are also possible, of course, but for simplicity we consider only the MTSI.) The most 
unstable mode has frequency o> ~ o) LH = < 0/(1 + w e 2 /O e 2 ) l/2 , with growth rate y ~ o) LH , and parallel wave number kn 
~ (m/M) l/2 k^. 

The salient property of the MTSI for our purposes is that it saturates by trapping ions in the perpendicular 
drift direction and electrons in the parallel direction; in the saturation process, the ions and electrons are heated to a 
fraction a 2 of the relative drift energy: 


Tj_i - T l|e - ot 2 MU ei 2 /2 . 


From simulations, McBride et al. (1972) find a — 0.5, with a wave energy density W at saturation of 
W/nMU ei 2 ~ a few percent. 
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Motivated by these results. Smith (1986b) postulates a self-scaling model in which the flank is assumed to be 
always at saturation (marginal instability) with respect to an instability such as the MTSI. The model characterizes 
the instability by two parameters a, (3, defined by 


V ii = a (cE x /B 0 ) ; U ei = (3 (cE x /B 0 ) , 


(4) 


where U ei is now the threshold drift speed and E x ~ Using equations (3) and (4), we find the self-similar 
scalings 


Rj 2 /C f 2 - a (3 ; T i ,/T eo = (a 3 (3) ,/2 (e<|) DL /T eo ) ; 



where X eo , T eo , and w eo are reference values of the electron Debye length, temperature, and plasma frequency. 

Owing to momentum conservation, there is a wave-modulated friction between the electrons and ions, 
which may be described by an anomalous collision frequency (Davidson and Krall, 1977) v* — ew LH , where e = 
W/nMUei 2 . Thus, the electron and ion fluids are acted on by volume forces F yi = -F ye , leading to an F x B drift 
velocity in the x-direction, i.e., opposite to E x (the coordinate system is defined by Fig. 3). This drift velocity is 
given by V xe = V xi = V x (n e ,<j) DL ), where 


3q5\1/4 


V x = e (a^) 


( m\ l/2 to LH / etfrpL^ 
M) ft, \T eo } 


1/2 


( 6 ) 


and V eo = (T eo /m) l/2 . Thus, above the DL, plasma is transported from the center of the arc to the flanks, concentrat- 
ing th e e° parallel current there (Fig. 3). Although V xe = V xi , there is a net current J x because above the region of 
strong E|| in the DL, we expect an extended region of small charge density p which sustains a weak parallel density 
field E|| driven by beam-plasma instabilities. Then the continuity equation is 


dJ z /dz — - dJ x /dx ~ - d(pV x )/dx 


Upon solving and integrating along the field line, we obtain (Smith, 1986b) 


J z (oo) _ J z (0) =» C(a,(3) (4) DL /ea 2 ) 5/8 


(7) 


where C(a,(3) is a constant and ft a is the perpendicular scale length of the arc (Fig. 3). The RHS of equation (7) is just 
the term f J x dz in equation (2). Assuming J z (0) « J z ( a0 ), equation (7) gives a scaling law 
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( 8 ) 


<W < (J z (°°)/c) 8/5 £ a 2 . 


For typical auroral-zone parameters equation (8) yields <J>dl ~ 10 kV for £ a ~ 1 km, in general accord with observa- 
tions. 


III. SIMULATION 


In the context of the above discussion, the flank may be modeled as a transmission line with local potential ct> 
= E x fi f , where £ f ~ <|> l/2 . Once the DL potential has reached a threshold value 4>* required to drive the MTSI, the 
transmission line is charged by the distributed (in z) perpendicular current J x = pV x of equation (7). We shall report 
elsewhere on simulations in which the transmission-line equations (Smith, 1986c) are solved simultaneously with a 
one-dimensional simulation of the DL evolution; this procedure provides the necessary self-consistent boundary 
condition on the current density Jdl( 0 at the simulation boundaries. In this paper, we replace the transmission line 
with a simple model circuit. 

If the flank were uniform between the DL (z = z* = 0) and the generator, the transmission line would appear 
to the DL as a pure impedance over the evolution time of the DL, with value Z-r = (Lr/C T ) I/2 , where Lr = £ f /4iTC 2 , 
C T = 4ir€j/£f. We model the impedance by the same form, with variable £f(<J>). We thus adopt the model circuit 
shown in Figure 4 where the diode symbol represents the DL and the variable resistor R(<J>) represents leakage 
current in the flank; this term is modeled by using the same form for J x as derived above, but over the perpen- 
dicular scale length £ f instead of £ a , by integrating J x ~ n l/2 (z) over the length z K (t) = ct/V'e i . The dielectric 
constant is defined by 


e l ~ 1 + (“eo 2 /^ 2 ) + (a> jo 2 /D, 2 ) [(l-ap) 2 + a 2 (l+7 e )] 


(9) 


where y e = T|, e /T j_,. 

The heuristic definition [equation (9)] is such that the total energy stored in the dielectric is € i E x 2 /8'ir: the 
first term in [] represents the reduced ion drift speed, while the second term accounts for ion and electron heating by 
saturation of the self-scaled modified two-stream instability. 

For the purpose of testing the scaling of 4>dl with £ a and a, we adopt the philosophy that owing to the 
separation between the perpendicular scale lengths £ a and £ f , the flank may be represented by these circuit elements 
while the DL will evolve in an essentially one-dimensional fashion in the central region of the arc. The DL evolution 
is simulated with the one-dimensional Vlasov code described by Smith (1982b), replacing the circuit used there by 
that of Figure 4. The boundary condition on the current density J DL (t) is then given by 

. / t x _ i s _ 4 ^dl(o r 1 1 1 

£ a £ a [Zr(A> DI ) + R(4> dl ) J 1 

where Is is the constant source current and 


Zt = [a/p] 1/4 (Y/e/ 72 ) (V eo /c) (£ a /p. 172 ) <|> DL ,/2 , 


( 10 ) 


( 11 ) 
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(12) 


R ■■ /V T eo\' /2 (Y£ a ) 2/3 <I>DL I/4 

\ 2 T J (a 7 p) ,/8 [ 1 -exp(-K 2 Z K (t))/2] 

|x = m/M, K 2 is defined in Smith (1986b) Y = a> eo /ft eo , V eo = 0 ) eo K eo , and all lengths and time scales are norma- 
lized to the nominal upstream quantities \ eo and a> eo , respectively; also 4 >dl->- e <J>DL/T eo . The term in [] in equation 
(10) replaces the physically derived I x of equation (8). This term is valid only after the threshold potential 4).+ has 
been attained, and so is turned on adiabatically for <j> > <|>*. Therefore, the circuit model does not accurately 
describe the initial dynamics in the linear instability phase of the evolution. In addition, the lumped circuit of Figure 
4 cannot represent the distributed nature of the flank charging, and so we cannot construct a circuit topology that 
allows for inductive fields. Therefore, we cannot model the acceleration of the inflowing (injected) distributions by 
inductive effects. 

For the parameters we use (see below), we estimate that the effect of neglecting inductive effects is small. As 
for the first limitation, the transient charging mechanism vanishes as d<j) DL /dt 0. Thus, we expect the model to be 
adequate for our present objective of studying time-asymptotic scalings. 

We show results for five runs. For all cases, the injected distributions are drifting Maxwellians with drift 
speeds in the simulation frame of U e = 2 V eo , Uj = -0.5 V eo . The forms of these distributions are held fixed (up to 
normalization). The threshold drift parameter (3 is held equal to 2, and M/m = 16. Holding C a /\ e = 60, we use 
values of a = 0.05, 0.02, 0.50. Fixing a = 0.50, we use S. a /\ e = 20, 40, 60. Initialization and other im- 
plementations are as described by Smith (1982b). 

Figure 5 shows the scaling of <|> DL with a, the fundamental parameter of the self-scaling marginal stability 
model of the MTSI discussed earlier. In the circuit equation (10), the principal effect of a is contained in the depen- 
ence of the impedance Z-r on the dielectric constant €^. In terms of the circuit equation (10), the DL scaling law 
equation (1) becomes in dimensionless notation 

. r. 4*dl ,/2 w ' 4 i 

*DL - GSolJ, , (13) 

where Zt(<J>) = Z 0 <J> I/2 , R(4>) = R 0 <J> I/4 . The RHS (13) has the form of a large factor GJ S £ DL 2 times a small factor 
[...], and the upper bound for <J> DL is obtained from setting [...] = 0. Because Z„ < < R 0 , equation (13) implies <j> DL 
~ Z G 2 . In Figure 5 we also plot the dependence of Z Q 2 on a, which agrees well with the plotted points. 

Figure 6 confirms the scaling <|> DL ~ £ a 2 found above. Again, this result is contained in equation (13) through 
the dependence <j) DL ~ Z c 2 (the results of Fig. 6 are all for a = 0.5, where Z 0 « R 0 ). 

Because the speed of light c is introduced in the impedance, the choice of V eo /c yields a physical scaling of 
velocities. Because Y = <a e jCl eo is a parameter, we obtain physical values of the length scales for an assumed value 
of either B c or n eo . All runs discussed here are for Teo = T io = 1 keV, typical of the plasma sheet population (note 
4 >du scales independently of T eo ). The scaled <|> DL is then given in kV as shown in the right-hand scale of Figure 6. 
Similarly, if we adopt a nominal value of B 0 = 0.05 G (f ce = 250 kHz) for the acceleration region, the top scale of 
Figure 6 gives arc thickness projected into the ionosphere of the order 1 km, which is the correct order of magnitude. 
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The power flows (I<|>) through the various parts of the circuit are shown for one case in Figure 7; for this case 
4 > dl = 42.5 kV. In the steady state, only 10 percent of the power flows goes through the DL, and about 90 percent 
goes into charging the flanks. The physical scale on the right shows the power dissipated per 1000 km extent of the 
arc in the E-W direction. Arcs are generally observed in systems of parallel bands, quasi-periodic in the N-S direc- 
tion; for the parameters of this example, each DA in such a system would dissipate about 10 1 1 W, compared with a 
typical substorm power of ~ 10 12 W. 

The scalings in a and C a have straightforward physical intepretations. The increase of <j> DL with a has two 
related aspects. First, the efficacy of the anomalous transport mechanism reported in Smith (1986b) increases with 
a, which is a measure of the strength of the MTSI. Second, in this self-scaling model the ratio Rj(<J>)/£f(<J>) ~ a 1/2 , so 
that as a increases, finite-Larmor radius effects lead to decreasing ion drift speed; hence, a higher ratio of electro- 
static to kinetic energy is stored in the flank “dielectric” for a given charge. The factor £ a in Z c [equation (6)] 
originates in the current balance (I s = J s £a)> arid because the flank is charged from the interior of the arc the charge 
available increases with £ a . The quadratic scaling <t> DL ~ 8 a 2 derives from the self-scaling of the transport model 
because £ f ~ <J> 1/2 . 


IV. DISCUSSION 


We have shown that basic considerations of DL evolution and stability require anomalous transport 
processes to divert the uniform upstream current to the flanks of a DA even after the parallel electric field has 
evolved to a steady state. The transport model we have discussed, albeit highly simplified, yields an estimate [equa- 
tion (8)] for the arc potential in general accord with observations. Other important consequences of the model are 
also in accord with satellite and rocket observations of DA’s and laboratory DL experiments. These include: (1) the 
density in a DA is substantially depleted relative to the ambient density (Benson and Calvert, 1979; Alport et al., 
1 986); and (2) concomitant with the depletion of the arc is that the current is diverted to the flanks, so that the highest 
current density is at the edges (Bruning, 1983; Burke, 1984). 

Besides the simple transport mechanism discussed here, there are many other mechanisms which are prob- 
ably important in DA’s. We are presently investigating models including ion-cyclotron modes. 

In terms of the simple circuit model, the potentials, perpendicular length scales, and power flows physically 
scale to correct orders of magnitude. For the nominal parameters we have chosen, the potential ranges from 5 to 42 
kV, while the length scales are consistent with the observational bound of —3 km on the latitudinal scale projected 
in the ionosphere (Boehm and Mozer, 1981). These quantities scale as n eo and ne 0 ' l/2 , respectively. We adopted B 0 
= 0.05 to correspond to the frequency of peak intensity of the auroral kilometric radiation, and chose w eo /fl eo as the 
marginal limit of strong magnetization (w e /l I e = R^/A.,,), which experiments reveal to be a requisite for strong DL 
formation with Earthward-directed Poynting flux (Smith, 1986a). (Recall that ne C is the ambient density before DL 
formation, not that of the evacuated arc.) Thus, our choice of co eo = fl eo is an upper bound; smaller values lead to 
lower potential, larger widths, and lower power flows. 

In future publications we shall report on refinements and extensions of the simulation concept, including a 
model in which the simplified circuit used here is replaced by the transmission-line equations. 
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Figure 1 . Schematic illustration of initial charging of the DL flanks by expulsion of charge from the localized, 
rapidly changing non-neutral region (shaded: stippled region p > 0, cross-hatched p < 0) where onset of current- 
driven instabilities occurs in the parallel sheet of a kinetic Alfven wave. 
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Figure 2. Schematic of the double layer flank spreading along B from the double layer toward the generator (here 
for illustration taken to be in the plasma sheet). 
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Figure 3. Left side: Schematic of the current diversion in the discrete arc. Right side: Definition of the scale lengths 
of the arc (£ a ) and flank (£ f ) and the associated scale factors which are used in the text to replace perpendicular 
derivatives. Also shown is a sketch of the inhomogeneous electric field E x , which produces relative drift between 
the electrons and ions owing to finite-Larmor-radius effects, and the definition of the xyz coordinate system. The 
magnetic field B 0 - -B c e z . 



Figure 4. Model circuit used to provide the current density boundary condition in one-dimensional DL simulation. 
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e<t>DL/T 


PROJECTED IONOSPHERIC 
DIMENSION* (km); B 0 = 0.4G 
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Figure 6. Scaling of 4> DL with the perpendicular scale width fi a of the arc model. The bottom and left scales are 
dimensionless. The right scale shows the potential in kV for assumed T eo = 1 keV; the top scale shows the arc 
dimension (2£ a ) projected into the ionosphere, for the given ambient parameters (before DL formation). 
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<|>dl (kV) 



RUN 8606 CIRCUIT POWERS 




Figure 7. Time history of power flows through the various circuit elements for run 8606 (e<t> DL /T eo = 42.5). TL - 
transmission line (flank) impedance; DL - double layer; AR - anomalous resistivity (leakage) in flank. The right- 
hand scale shows the physically scaled power for an arc extended 1000 km in the E-W direction. 
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